Complex lithium ion dynamics in simulated LiPC>3 glass studied by means of 

multi-time correlation functions 
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Molecular dynamics simulations are performed to study the lithium jumps in L1PO3 glass. In 
particular, we calculate higher-order correlation functions that probe the positions of single lithium 
ions at several times. Three-time correlation functions show that the non-exponential relaxation of 
the lithium ions results from both correlated back-and-forth jumps and the existence of dynamical 
heterogeneities, i.e., the presence of a broad distribution of jump rates. A quantitative analysis yields 
that the contribution of the dynamical heterogeneities to the non-exponential depopulation of the 
lithium sites increases upon cooling. Further, correlated back-and-forth jumps between neighboring 
sites are observed for the fast ions of the distribution, but not for the slow ions and, hence, the back- 
jump probability depends on the dynamical state. Four-time correlation functions indicate that an 
exchange between fast and slow ions takes place on the timescale of the jumps themselves, i.e., the 
dynamical heterogeneities are short-lived. Hence, sites featuring fast and slow lithium dynamics, 
respectively, are intimately mixed. In addition, a backward correlation beyond the first neighbor 
shell for highly mobile ions and the presence of long-range dynamical heterogeneities suggest that 
fast ion migration occurs along preferential pathways in the glassy matrix. In the melt, we find no 
evidence for correlated back-and-forth motions and dynamical heterogeneities on the length scale of 
the next-neighbor distance. 

PACS numbers: 66.30.Dn 



I. INTRODUCTION 

In view of numerous technological applications and a 
large variety of phenomena, recent scientific work repeat- 
edly focused on glassy ion conductors. On the one hand, 
it is well known that the macroscopic charge transport in 
these materials results from the diffusion of mobile ions 
in a basically rigid glassy matrix, on the other, a detailed 
microscopic picture of this dynamical process is still lack- 
ing. One key feature of the dynamics in glassy ion con- 
ductors is the non-exponential relaxation of the mobile 
ions, indicating the complexity of the motion. For in- 
stance, broad loss peaks are found in dielectrical and me- 
chanical relaxation experimentsiSii^. Moreover, when 
the dynamics of the mobile ions are studied by means of 
advanced nuclear magnetic resonance (NMR) techniques, 
non-exponential correlation functions are observe d-'?'^'" . 
In general, the origin of non-exponential relaxation can 
be twofold^. In the homogeneous scenario, the dynam- 
ics of all particles are characterized by the same relax- 
ation function that is intrinsically non-exponential due 
to a specific mechanism of the motion, e.g., correlated 
back-and-forth jumps. Contrarily, in the heterogeneous 
scenario, all particles relax exponentially, i.e., there are 
no correlated back-and-forth jumps, but a distribution of 
correlation times G(lgr) exists. In our case of a complex 
motion, one may expect that both homogeneous and het- 
erogeneous dynamics contribute to the non-exponential 
relaxation. 

A strong frequency dependence of the electri- 
cal conductivity cr(u>) of ion conductors indicates 
that back-and-forth motions occur during ionic 
diffusioniiii2iiiiiiSiiSiii. However, to study a complex 



motion in detail, it is not sufficient to probe the dynamics 
at two times, but one has to resort to multi-time corre- 
lation functions^. For example, three-time correlation 
functions can be used to quantify the extent to which 
the non-exponential relaxation results from the homo- 
geneous and the heterogeneous scenario, respectively^. 
Four-time correlation functions allow one to measure the 
lifetime of dynamical heterogeneities, i.e., the timescale 
of exchange processes between fast and slow ions. Two 
methods have recently proven well suited to record 
multi-time correlation functions for glassy ion conduc- 
tors. While ion dynamics on the ps — ns timescale can 
be studied in molecular dynamics (MD) simulations^, 
multi-dimensional NMR techniques probe jumps on a 
timescale of ms — s. Specifically, multi-dimensional 109 Ag 
NMR experiments showed that a broad rate distribution 
governs the jumps of silver ions in silver phosphate 
based glasses! Hence, both correlated back-and-forth 
jumps and dynamical heterogeneities indeed contribute 
to the non-exponentiality. However, to the best of our 
knowledge the relevance of both contributions has not 
yet been quantified for glassy ion conductors. 

MD simulations are an ideal tool to investigate dy- 
namical processes on a microscopic level. Prior work 
on the origin of the non-exponential relaxation in glassy 
ion conductors has focused on alkali silicate glasses. For 
lithium silicate glasses, it was reported that dynamical 
heterogeneities contribute to the non-exponential relax- 
ation of the lithium ions where the lifetime of these het- 
erogeneities is limitedi2*2&. Further, back-and-forth mo- 
tions on various length scales were observed, including a 
backward correlation beyond the first neighbor shell at 
low temperatures^. For sodium silicate glasses, no evi- 
dence for back-and-forth jumps between adjacent sodium 



2 



sites was found at higher temperatures^. Moreover, sev- 
eral workers demonstrated that, though there is no micro- 
segregation, the sodium ions follow preferential pathways 
in the glassy matrix 22 '—. The sodium ions inside these 
channels show a higher mobility than the ones outside 
and, hence, dynamical heterogeneities exist 24 . In addi- 
tion, studying examples of ionic trajectories the mech- 
anism of alkali ion migration in silicate glasses was de- 
scribed as "vacancy-like process"—. During a correlated 
motion of two cations, an alkali site is left by one ion an 
occupied by another. 

In the present work, we use MD simulations to ana- 
lyze the complex mechanism of ion dynamics in an alkali 
phosphate glass. In this way, we intend to complement 
the prior studies on silicate glasses and to compare the 
nature of ion diffusion in different amorphous materials. 
Moreover, based on an analysis of novel higher-order cor- 
relation functions, new information about this dynami- 
cal process is revealed. In detail, the dynamics of non- 
crystalline LiPC>3 is investigated in a temperature range 
where the motions of the various atomic species strongly 
decouple upon cooling. While the lithium ionic sub- 
system can be equilibrated at all chosen temperatures, 
the diffusion of the phosphorus and the oxygen parti- 
cles freezes in on the ns-timescale of the simulation. In 
other words, the dynamics in the melt and in the glass 
are studied. Our main goal is to identify the tempera- 
ture dependent origin of the non-exponential relaxation 
of the lithium ions. For this purpose, we calculate multi- 
time correlation functions that link the positions of sin- 
gle ions at subsequent points in time. In this way, we 
quantify the homogeneous and the heterogeneous contri- 
butions to the non-exponential relaxation and measure 
the lifetime of the dynamical heterogeneities. Since the 
multi-time correlation functions observed in NMR exper- 
iments on solid ion conductors and some of the quantities 
computed in the present study have a similar information 
content, a future comparison of the respective results is 
very promising. 



II. DETAILS OF THE SIMULATION 

The potential used to describe the interaction of the 
ions in LIPO3 can be written as the sum of a Coulomb 
and a Born-Mayer-Huggins pair potential 

<$> a p{r) = q -^- + A a0 exp {-r/p) (1) 

where r is the distance between two ions of type a and 
type /3, respectively. The potential parameters, except 
for Ano, are adopted from the work of Karthikeyan et 
al— . Here, Ano is reduced to obtain a more realistic 
interatomic distance rno, cf. below. In detail, we use 
effective charges qn = 0.6, qp — 3.0 and q = —1-2 as 
well as p = 0.29A. The parameters A a p are listed in Tab. 
1. Most simulations are performed in the NEV ensemble 



TABLE I: Interaction parameters A a f3, cf. Eq. 1, together 
with the interatomic distances r a p as obtained from the 
present simulations and from experimental work—*—, respec- 
tively. 





A Q/3 [eV] 


r Q/ 3 [A] 
simulation 


T a /3 [A] 

experiment 


LiLi 


167.47 


2.74 




LiP 


158.09 


3.31 




LiO 


300.00 


2.03 


2.02 


PP 


148.93 


3.33 


3.01 


PO 


694.66 


1.51/1.66 


1.52-1.56 


OO 


2644.8 


2.58 


2.52 



for N — 800 particles. To take into account the ther- 
mal expansion of phosphate glasses, the density is hxed 
at p = 2.15 g/cm 3 as compared to p — 2.25 g/cm 3 ob- 
served in experiments at room temperature^. This ap- 
proximation for the density at higher temperatures is ob- 
tained based on the linear expansion coefficient of LLPO3 
glass^. In addition, we consider systems with N — 200 
and N = 400 to check for finite size effects. Moreover, 
some simulations are performed at "constant pressure" 
to enable a comparison with experimentally determined 
activation energies. Strictly speaking, we continue doing 
runs in the NEV ensemble, but the size of the simulation 
box is adjusted so that a pressure p = 7.0 ± O.lGPa is 
obtained at all temperatures. This procedure results in 
a decrease of the density from p — 2.21 g/cm 3 at 744K 
to p = 1.95 g/cm 3 at 3000K. The equations of motion 
are integrated using the velocity form of the Verlet algo- 
rithm with a time step of 2fs. Further, periodic boundary 
conditions are applied and the Coulombic forces are cal- 
culated via Ewald summation. In doing so, consistent re- 
sults are obtained when computing the trajectories with 
the programme MOLDY22i and an own MD simulation 
programme, respectively. While the whole system can 
be equilibrated at T > 1006K, the phosphate glass ma- 
trix is basically frozen on the timescale of the simulation 
at lower temperatures. In the latter cases, we pay atten- 
tion that the lithium ionic subsystem is still in equilib- 
rium. For example, the sample is equilibrated for 20ns 
at T = 592K before recording the data. 

Karthikeyan et al— demonstrated that the structure 
of LiPC>3 glass is well reproduced when using the two- 
body potential specified in Eq. ^ We confirmed this 
result by analyzing the structure of the glass obtained 
after a quench from T = 592K to T = 300K. Tab. 1 shows 
that the interatomic distances r a p are in good agreement 
with the corresponding experimental values 2 ^— Typical 
of phosphate glasses^!, the P-0 distance for tangling oxy- 
gens (rpTo — 1.51 A) is shorter than for bridging oxygens 
{fpBO = 1.66 A). Further, the reduction of A^o, as com- 
pared to the value used by Karthikeyan et al^, improves 
the agreement with the actual distance r^o^Li^S. Con- 
sistent with experimental findings for LiP03 glass^S 2 ^, 
the bond-angle distributions and the coordination num- 
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FIG. 1: (a): MSD of the lithium ions, r|,(t), in non- 
crystalline LiPOs at the indicated temperatures. Dashed line: 
r 2 Li{t) = GD Li t with D u =8.31 -10" 7 cm 2 /s. (b): Temperature 
dependence of the diffusion constants Dli, Do and Dp; solid 
symbols: TV = 800, p = 2.15g/cm 3 , open symbols: iV = 800, 
p = 7.0GPa, vertical stripes: TV = 400, p = 2.15g/cm 3 , horizon- 
tal stripes: TV = 200, p = 2.15g/cm 3 , lines: Arrhenius fits to 
the data for TV = 800. 
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FIG. 2: Incoherent intermediate scattering function F a {q m = 
2.28A ,t) for the lithium ions in non-crystalline LiPOs at 
the indicated temperatures (TV = 800, p — 2.15g/cm 3 ). Inset: 
Temperature dependence of the mean correlation time tfs', 
solid symbols: TV = 800, vertical stripes TV = 400, horizontal 
stripes: TV = 200, line: Fit to an Arrhenius law for TV = 800. 



bers indicate that the simulated glass consists of well de- 
fined phosphate tetrahedra that are connected by two 
of their corners to form long chains and/or rings. Con- 
cerning the intermediate range order, some differences 
between simulated and actual LiPOs glass may exist. In 
agreement with results by Karthikeyan et ali^i, we find 
a mean P-O-P bond angle that is smaller than the one 
observed in experiments on vitreous P 2 0g2i. To fix this 
mean angle at the experimental value for vitreous P2O5, 
Liang et al^ used a three-body potential in simulations 
of L1PO3 glass. We refrain from doing so because, first, 
the P-O-P bond angle in meta-phosphate glasses has not 
yet been determined experimentally and, second, three- 
body interactions distinctly slow down the simulation. 
For our analysis, it is important to perform simulations 
at sufficiently low temperatures, because not until then 
lithium dynamics show several characteristic features, 
e.g., back-and-forth jumps between adjacent site o 19 ' 35 . 
Therefore, we keep the potential as simple as possible so 
as to be able to equilibrate the system at these temper- 
atures. 

III. RESULTS 

To study the dynamics of non-crystalline L1PO3, we 
first display the mean square displacement (MSD) of the 
lithium ions 

rlS)= {[r j {t)f)= ({\ri(t + t)-r>(t )\} 2 ) (2) 

in Fig. |IJa). Throughout this article, the brackets (. . .) 
denote the average over the lithium ionic subsystem, r J (t) 
is the displacement of lithium ion j during a time inter- 
val t and f J (to) specifies its position at the time to. At 
temperatures T<2040K, r^^t) exhibits three character- 
istic time regimes: ballistic motion r 2 Li ott 2 at short times 



t < O.lps, sublinear diffusion due to back-and-forth mo- 
tions of the lithium ions at intermediate times and linear 
diffusion r 2 Li rxt at long times. Upon cooling the diffusion 
of the lithium ions slows down and the time window of 
sublinear diffusion is extended. A qualitatively similar 
behavior is found for the MSD of the oxygen and the 
phosphorus ions, ^(rj) and r 2 P {t) (not shown). 

The diffusion constants Duo.p extracted from the 
long-time behavior of r 2 Li p {t) are compiled in Fig. 
[Jb). Evidently, the lithium diffusion decouples from 
the oxygen and the phosphorus dynamics. While Dn 
follows an Arrhenius law with activation energy Ed = 
0.40eV, the curves lg Do,p(l/T) show a downward bend- 
ing. In particular, on the lOns-timescale of our simula- 
tion, a basically rigid phosphate glass matrix exists be- 
low T w 1000K. Therefore, we refer to T g = 1000K as 
the computer glass transition of LiPC>3. For comparison, 
the diffusion constants resulting from the calculations for 
TV = 200 and TV = 400 and from the constant pressure sim- 
ulations are also included in Fig. 1(b). We find no indi- 
cation that the diffusion constants depend on the system 
size. On the other hand, the temperature dependence 
of Du in the constant pressure runs {Ed = 0.47eV) is 
stronger than in the constant volume simulations. 

Next, the incoherent intermediate scattering function 
F s {q,t) for the lithium ions is considered: 

F s {q,t) = < cos [q{r 3 {t +t) - r> {t ))}> . (3) 

In our case of an isotropic sample, this function only de- 
pends on the absolute value, q, of the wave vector. To ob- 
tain information about dynamical processes on the length 
scale of the Li-Li interatomic distance, q = q m = 2Tr /r^iu 
is used. Fig. |21 shows F s {q m , t) for the NEV ensem- 
ble (TV = 800). We see temperature dependent, non- 
exponential decays. These findings can be quantified by 
fitting the data to a Kohlrausch- Williams- Watts (KWW) 
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FIG. 3: Representative trajectories of the lithium ions in non- 
crystalline LiPOa. Left hand side: Trajectories at T = 2040K 
during a time interval At = 3ps. Right hand side: Trajectories 
at T = 592K during a time interval At = 100ps. 



function^, A exp[— (t/r) 13 ]. Then, a stretching param- 
eter (3 rs 0.42 characterizes the non-exponentiality at 
T < T g . Above T g , an increase to (3 « 0.54 is ob- 
served, however, the absence of a short-time plateau of 
F a (q m ,t) leads to larger error bars in this temperature 
range. In the inset of Fig. we display the mean time 
constants Tps calculated from the fit parameters accord- 
ing to Tps — ( r //?)r(l//3) where T(x) is the T-function. 
Though there are some deviations, the temperature de- 
pendence of tfs can still satisfactorily be described by 
an Arrhenius law with activation energy Eps — 0.62eV. 
In the same way, Eps = 0.67 ± 0.03eV is obtained from 
the constant pressure simulations. Thus, in both cases, 
Eps is significantly larger than Ep>. In addition, Fig. 2 
shows that the time constants Tps do not depend on N 
and, hence, finite size effects are again absent. In the re- 
mainder of this article, we analyze the data from constant 
volume simulations for TV = 800. 

To study the mechanism of lithium dynamics, we first 
look at some representative trajectories. Fig. [3] displays 
projections on the xy-plane for trajectories during com- 
parable time intervals at T = 2040K and T = 592K. While 
the mechanism of the motion resembles liquid-like dif- 
fusion at the high temperature, it is rather complex in 
the glass. Specifically, the following features are obvi- 
ous for T = 592K: (i) The motion of most lithium ions 
can be decomposed into local vibrations and jumps be- 
tween adjacent sites, (ii) The lithium ions show differ- 
ent mobilities, (iii) Back-and- forth dynamics involve dis- 
tinct lithium sites, (iv) The number of correlated back- 
and-forth jumps varies. While numerous back-and-forth 
jumps are observed for ion 4, only a few are found for 
ion 3. (v) Several ions (6-10) follow the same preferen- 
tial pathway. Consistent with the results of simulations 
on alkali silicate glasses^i^^ 3 ^ 9 ., all these features elu- 
cidate the complexity and the diversity of ionic diffusion 
in glasses. 

A statistical analysis of the mechanism of the motion is 
possible based on the self part of the van Hove correlation 



FIG. 4: Space- and time dependence of the self part of the van 
Hove correlation function G„ (r, t) for non- crystalline LiPOa at 
the indicated temperatures. 

function for the lithium ions. It is given bj^S 

G s (r,t) = (S(r-r j (t))) (4) 

where S(x) denotes the ^-function. G s (r,t) measures the 
probability that a lithium ion moves a distance r in a time 
interval t. Fig. Q] shows G s (r,t) for three characteristic 
temperatures. For T = 2040K, we observe a single maxi- 
mum that shifts to larger values of r when t is increased. 
Though the findings resemble expectations for liquid-like 
diffusion, an asymmetry of the curves for t ks lps indi- 
cates deviations. For T<T g , G s (r, t) exhibits oscillations 
that become more pronounced upon cooling. Specifically, 
there are minima at rssl.7A and r«4.2A together with 
a maximum at r «2.7A«r£i£i. Thus, the glassy network 
provides well defined lithium sites that are separated by 
energy barriers E > ksT so that, to a good approxi- 
mation, the lithium ionic motion can be described as a 
sequence of hopping processes, i.e., the timescale it takes 
to cross the barriers is much shorter than the timescale 
the ions fluctuate about the sites. In particular, it is 
possible to distinguish between ions that, after a time t, 
(i) still reside at the same site < 1.7A), (ii) have 

jumped to a next-neighbor site (1.7A< r J (<) <4.2A) and 
(iii) have moved beyond an adjacent site (r J (i) >4.2A). 
These findings for G s (r, t) agree with both the appear- 
ance of the trajectories in Fig. [3] and the results of pre- 
vious simulationa 19 . 35 ! 39 . 41 ! 42 . 43 ! 44 . 

Next, we define the two-time correlation function 
S 2 {t) = (s|(t)) where s|(<) = 1 for < l - 7 ^ and 

s i W = otherwise. In our case of hopping motion, S2 (t) 
measures the probability that a lithium ion is still at the 
initial site after a time t. Thus, it directly reflects the de- 
population of the initially occupied lithium sites. In Fig. 
[SI we see that ^(t) decays non-exponcntially. This non- 
exponentiality can result from correlated back-and-forth 
jumps between neighboring sites and/or from dynami- 
cal heterogeneities, i.e., the presence of a distribution of 
jump rates. Based on an appropriate three-time corre- 
lation function S3, both contributions will be quantified 
later in this article. The average time it takes for the 
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FIG. 5: Two-time correlation function Sity) at the indicated 
temperatures. 



lithium ions to successfully escape from their sites, i.e., 
the timescale of the depopulation of the sites, can be 
characterized by the 1/e-decay time of S 2 (t), T S2- Fig- 
El shows that the temperature dependence of tsi is well 
described by an Arrhenius law with activation energy 
S 5 = 0.44eV. 

In the remainder of this article, the origin of the non- 
exponential relaxation of the lithium ions in simulated 
LiP03 glass is studied in detail by analyzing various 
multi-time correlation functions. All these functions cor- 
relate the positions of single lithium ions at subsequent 
times t\<t 2 < .... To characterize the motion of ion j be- 
tween two times t a <tp, we use the notation r£ =f J (t a ), 
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'a/3 = ' /3 ~'a> ' Q/ 3 = 1' a/31 t ap= t (3 

To show the existence of dynamical heterogeneities and 
to study their lifetime, we define a "four-time van Hove 
correlation function" 

Gi{r,ti2=t s ,t23,hi = ts) = {a(r 12 ) $(r - r^)). (5) 

In this equation, the term 0(7-12) simply means that three 
subensembles are selected with respect to the displace- 
ment during the time interval t±2 = t s : SEo consists of 
ions that are still at the initial site (r^ (tn) < 1-7A), 
SEi contains the ones that have jumped to an adja- 
cent site (1.7A < r J '(ii 2 ) < 4.2A) and the ions of SE 2 
show r J (ti2) > 4.2A. Thus, Gi(r, t s ,t23,t s ) probes the 
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FIG. 6: Various time constants characterizing lithium ion dy- 
namics in non-crystalline LiPOs, see text for details. 




FIG. 7: Four-time van Hove correlation function Gi(r, t s ,t, t 3 ) 
for the lithium ions in L1PO3 glass at T = 713K. A filter time 
£ 3 = 204.8ps was used to select three subensembles of lithium 
ions SEq, SEi and SE2, see text for details. 



van Hove correlation function G s (r,t s ) of the respective 
subensembles starting at a time t 23 after their selection. 
The results for t s = 204. 8ps and various t 23 at T = 713K 
are displayed in Fig. El For ^23 = 0, the three subensem- 
bles show clearly different G4. In particular, the first 
peak is highest for SEo, indicating that ions that have not 
escaped from their sites during t\ 2 move also less then the 
average particle in the time interval £34. In other words, 
slow and fast ions are distinguishable and, hence, dynam- 
ical heterogeneities exisfeiS. An extension of t%3, i.e., of 
the delay between selection and detection, allows us to 
study the lifetime of the dynamical heterogeneities. In 
Fig. [7J we see that the curves for the different subensem- 
bles approach each other when £23 is increased until they 
are nearly identical after a time £23 = 4.6ns. At this time 
t = 4.6ns, the decay of S 2 (t) is just complete, cf. Fig. 
Thus, the information about the initial dynamical states 
is lost when all ions have successfully escaped from their 
sites, i.e., the dynamical heterogeneities are short-lived. 
A more detailed analysis of the lifetime of the dynamical 
heterogeneities will be presented below. 

For a study of back-and-forth jumps, we focus on 
lithium ions that have moved to a neighboring site during 
ti 2 and map their motion in the subsequent time interval 
/;23- This analysis is possible based on the "three-time 
van Hove correlation function" 



G 3 (r,t 12 ,t 23 ) = (d(r{ 2 )6(r-r{ 3 )) 



(6) 



where d{r{ 2 ) = 1 for ljA < r{ 2 < 4.2A and d{r{ 2 ) = 
otherwise. G 3 (r, tn, £23) measures G s (r, i = ti2 + £23) ex- 
clusively for ions that occupy a next-neighbor site at the 
time t 2 . G 3 (r, t\%, t 23 = 0) is different from zero only 
in the range 1.7A < r < 4.2A. For t 23 > 0, backward 
jumps towards the initial positions r/ yield the inten- 
sity lb at r < 1.7A, whereas forward jumps contribute 
to the intensity If at r > 4.2A. Hence, when multiple 
jumps are negligible, i.e., for £23 < £i2> the ratio lb/ If 
provides information whether the subsequent motion is 
backward or forward correlated. In Fig. |S] we display 
G 3 (r,h 2 ,t 23 ) at T = 592K for t 12 = 12.8ps, 819. 2ps and 
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Gz{r,tv2,tiz) for the lithium ions in L1PO3 glass at T = 592K r i2^ 



various *23- Nearly independent of £23 < £12, integration 
yields h/If ~ 15 for t%2 = 12.8ps and h/If ~ 0.3 for 
ti2 = 819. 2ps. While the former ratio indicates a sig- 
nificant backward correlation, the latter is comparable 
to h/If ~ 1/3 as expected for a random- walk based on 
a mean coordination number z^n^i w 4. For an inter- 
pretation of these results, it is useful to recall the exis- 
tence of dynamical heterogeneities and to reinspect Fig. 
0J It becomes clear that the population at the neighbor- 
ing site is still small after a time t\2 — 12.8ps and, hence, 
it mainly results from the fast ions. Vice versa, ions with 
a medium or a low mobility occupy an adjacent after a 
time tu — 819. 2ps. Hence, the curves G 3 {r, i i2 , £23) for 
t\2 = 12.8ps and t\2 = 819.2ps reflect the behavior of fast 
and non-fast ions, respectively. We conclude that a sig- 
nificantly enhanced back-jump probability exists for the 
fast ions, but not for the slow ions. 

One may argue that large-amplitude vibrations can 
also lead to h/If 3> 1 for short t 12 - However, one 
expects that displacements due to vibrations are lim- 
ited by r\ 2 < Thihi and, hence, a significant contribu- 
tion of this type of motion would result in an asym- 
metry of G 3 (r, ti2, £23 — *-0). Therefore, the nearly sym- 
metric shape of G 3 (r, £12, £23 < 0.4ps) shows that large- 
amplitude vibrations do not affect our conclusion. The 
finding that the back-jump probability depends on the 
dynamical state is backed up by preliminary results of 
a more detailed analysis. There, we explicitly identify 
the lithium sites and calculate the back-jump probabil- 
ity pb as a function of the waiting time at the initial 
site, t w . Consistent with results for LiSi0 2 45 , we find 
that Pb(t w ) significantly decreases with increasing wait- 
ing time. Thus, a high backward correlation exists exclu- 
sively for ions that show a high jump rate at the initial 
site. 

The length scales of the dynamical heterogeneities and 
of the correlated back-and-forth motions can be studied 
based on the conditional probability functions p (2^23 1 ^12) 
and p (2/23 1 ^12)- These functions measure the probabil- 
ity to find specific values x 3 ^ = X23 and y 23 = t/23, re- 
spectively, provided the particle has moved a distance 
r i2 = r i2 m the first time interval t\2- Here, x J 23 is the 



FIG. 9: First moment £23 (^12) of the conditional probability 
function p (X23 \ r 12) for LiPOs glass at T — 592K. The time 
intervals ti2 and t23 are indicated. 

projection of r 2 3 on the direction of the motion dur- 
ing tu, i.e., x 3 23 = r 2 \ ■ rl 2 /r[ 2 and y J 23 is the projec- 
tion of f 23 on an arbitrary direction perpendicular to 
rf 2 . Motivated by the outcome of prior work 19-46,47 , we 
focus on the first moment X23 as well as on the sec- 
ond moments a x (ri 2 ) =< [x J 23 -X23{ri2)] 2 (r{ 2 ) > and 

o-j,(ri 2 ) = < [2/23-1/23 ( r 12)] 2 (r{ 2 ) >, rather than analyzing 
the full probability functions. 

For the first moment, x 2 3(ri 2 ) = will result if 
the directions of the motions during t±2 and £23 are 
uncorrelated 4 ^ 7 .. In contrast, if the subsequent motion 
for a given r\ 2 is backward (forward) correlated, a neg- 
ative (positive) value of 2:23 (ri 2 ) will be observed. In 
particular, it has been readily shown that x 2 3(ri 2 ) = 
— (l/2)ri2 is found for stochastic dynamics in a harmonic 
potential 19 . Fig. shows £23(7*12) for T = 592K. For all 
t\2 and ^23, 5T23(ri2) follows the curve — (l/2)n 2 up to 
7*12 = 1.0 — 1.5A, indicating that the lithium sites can be 
described by harmonic potentials on this length scale. At 
larger n 2 , ^23 (^12) depends on both t\2 and £23. For a 
given ii 2 , X23(ri2) is more negative the lager the value of 
£23 until a saturation is reached at £23 ~ lOOOps. This can 
be understood as follows: No significant motion occurs 
during short time intervals £23 so that x 3 23 ~ 0. How- 
ever, when i 2 3 is extended, more and more lithium ions 
show a distinct displacement until the average behavior 
of all ions with a given r\2 is observed for long £23- To 
study this average behavior we further discuss the re- 
sults for a long £23 = lOOOps. In this case, x~23(?\2) is 
basically constant at ri 2 > lA for £12 = lOOOps, whereas, 
for t\2 < lOOps, a decrease over the whole accessible ri 2 
range indicates that the backward correlation increases 
even beyond the interatomic distance r^Li. A compari- 
son with Fig. 0] shows that, for ti2 < lOOps, the findings 
in the range run < r\2 < 8A, say, reflect the behavior 
of fast ions, while an ion with an average mobility shows 
such displacements during a time interval t±2 — lOOOps. 
Thus, the decrease of x 2 3(ri 2 ) at r\2 > ruLi for short 
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FIG. 10: Second moments a x (ri2) (solid symbols) and a y (r 12) 
(open symbols) of the the conditional probability functions 
p{x2z\i'i2) and p (j/23 1 fii), respectively. Data for non- 
crystalline LiP0 3 at T = 2040K and T = 592K are compared. 
The time intervals t — ti2 = t23 are indicated. 



ti2 implies that the correlated back-and-forth motions 
of fast ions may involve not only one, but also several 
sites. On the other hand, X23 ~ const, at ri 2 > lA for 
long ti2 again shows that the non-fast ions perform local 
rattling motions, but no correlated back-and-forth jumps 
between adjacent sites. We add that X23 ~ — 0.24ri 2 at 
ri2 < 2k and X23 = const, at r i2 > 2A are found for 
LiP0 3 melt at T = 2040K (t 12 < lps, t 23 > lps). Qualita- 
tively, these findings resemble the outcome of MD sim- 
ulations on supercooled liquids^^I. There, the results 
were interpreted as a signature of the cage effect describ- 
ing that, though all particles are mobile, every particle is 
temporarily captured in a cage built up by its neighbors. 
Then, back-and-forth motions take place within the cage, 
but not when the particles have left their cages. 

Now, we turn to the second moments. If all parti- 
cles have equal mobility, the displacements r±2 and r 2 3 
are uncorrelated so that u x , y (r\2) = const.. In contrast, 
if dynamical heterogeneities exist, on average, slow and 
fast ions show small and large ri 2 , respectively. Hence, 
ions with small ri 2 move less in i 2 3 than those with large 
r\2 and o- Xt y(r\ 2 ) increases. Fig. I1UI shows a x (r\2) and 
o~ y {r 12) for various time intervals t = t\ 2 = t 2 3 at two 
temperatures. For T = 2040K, the second moments are 
equal and independent of r i2 for t > lps, indicating the 
absence of dynamical heterogeneities related to lithium 
diffusion. In contrast, <7 x (ri 2 ) and cr y (ri 2 ) for t < lps 
increase, though, with a different slope and, thus there 
are highly anisotropic dynamical heterogeneities on short 
timescales. Since rn(t < lps) < rnLi> cf. Fig. 1, we at- 
tribute the latter dynamical heterogeneities to the rat- 
tling motions within the local cages. For T = 592K, a 
strong r\2 dependence of a XiV is observed for all chosen t. 
This finding again implies the existence of dynamical het- 
erogeneities attributed to lithium diffusion in the glass. 



Since back jumps during t 2 3 result in x 23 7^0 and y 2 3 «0, 
i.e., the component perpendicular to the direction of the 
forward jump is small, the dynamical heterogeneities can 
be most directly studied based on o~ y (r\2)- The ratio 
o~ y (rLiLi)/o- y (0)^>l shows that the heterogeneity results 
for the most part from different jump rates at neighbor- 
ing lithium sites. However, additional contributions on 
a length scale of several Li-Li interatomic distances are 
indicated by a further increase of <J y {r 12) up to ri 2 « 13A 
for t > lOOOps. Considering also the appearance of the 
lithium trajectories in Fig. [21 we suggest that preferred 
pathways of ion migration connect several lithium sites. 

The oscillatory behavior of ir 2 3(ri 2 ) and o~ x , y {r 12) al- 
lows us to study the energy landscape governing the 
lithium dynamics. For X23{t\2) 1 one expects local min- 
ima (maxima) near positions where the slope of the po- 
tential shows maxima (minima) 1 ^. In Fig. [§1 we see that 
£23(^12) shows a minimum at ryi ~ l.lA and a maximum 
at Tyi ~ 2.3A. Hence, there should be an energy barrier 
at an average distance r w 1.7A. At r 12 « 1.7A, the curves 
@x. y (7*12 ) exhibit local maxima. This is reasonable, be- 
cause ions residing at a saddle after t i2 should possess a 
comparatively high mobility in the subsequent time in- 
terval. A closer inspection of the data reveals that the 
maxima of u x . v {r 12) shift to smaller ri 2 when t is ex- 
tended. As will be discussed later in this article, this 
finding is consistent with the assumption that the Li-Li 
interaction depends on the waiting time at a site. 

To quantify the homogeneous and the heterogeneous 
contributions to the non-exponential relaxation of the 
lithium ions we define the three-time correlation function 

S3(t 12 =t/2, t 23 =t/2) = <4(*i2)4(*2 3 )> ■ (7) 

For purely homogeneous and purely heterogeneous re- 
laxation, such three-time correlation functions can be 
expressed by the corresponding two-time correlation 
functiono 18 i 46 i 48 . 49 . While S 3 (t/2,t/2) = [S 2 {t/2)} 2 is 
valid in the homogeneous scenario, S3(t/2,t/2) = «S 2 (f) 
holds in the heterogeneous scenario, cf. Appendix A. In 
Fig. ^2 S f 3 (i/2, i/2) is shown for three temperatures. For 
T = 1287K, S 3 (t/2,t/2)w[S2(t/2)} 2 implies that the re- 
laxation is basically homogeneous. However, above T g , 
hopping dynamics is not a good approximation for the 
mechanism of the motion and the phosphate matrix dy- 
namics modifies the lithium sites. Therefore, our findings 
do not indicate back-and-forth jumps between adjacent 
lithium sites, but rather they reflect back-and-forth mo- 
tions that are affected by the matrix dynamics and take 
place on various length scales. At T < T g , S , 3 (i/2,</2) 
clearly deviates from [^(t/2)] 2 and, hence, dynami- 
cal heterogeneities become increasingly important upon 
cooling. Their relevance can be further analyzed based 
on the ratio 



H(t) = 



S3(t/2,t/2)~{S2(t/2)}< 



S 2 (t)-[S 2 (t/2)}2 
It is easily seen that H(t) = 1 and H{t) 



(8) 

hold for 
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FIG. 11: Three-time correlation functions Ssit/2, f/2) for 
non-crystalline LiP03 at the indicated temperatures. The 
results are compared to the expectations for purely heteroge- 
neous dynamics, 52 (t), and for purely homogeneous dynam- 
ics, [S 2 (t/2)] 2 . 
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FIG. 12: H(t) for non-crystalline L1PO3 at the indicated tem- 
peratures. The data were calculated from the results in Fig. 
11 using Eq.[H] 



purely heterogeneous and purely homogeneous dynam- 
ics, respectively. In Fig. ^1 H{t) is displayed for the 
studied temperatures. Evidently, the non-exponentiality 
at low temperatures results to a large extent from dy- 
namical heterogeneities, i.e., there is a broad distribu- 
tion of jump rates. For T< 1006-ftf, the curves H(t) show 
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FIG. 13: Four-time correlation functions S4,(t 3 ,t,t a ) for non- 
crystalline LiP03 (points) and fits to a KWW function 
(dashed lines), (a): Data at T = 713K for the indicated 
filter times t 3 . (b): Temperature dependence of the curves 
for comparable filter times, i.e., S 2 (t s ) = 0.48 ± 0.02 (592K: 
t s = 145.0ps, 713K: t s = 40.0ps, 1006K: t s = 4.0ps, 1287K: 
i s = 1.2ps). 



a maximum at a time th that increases upon cooling. 
Inspecting the temperature dependence of th in Fig. |HJ 
th{T)~ts2(T) is found. Thus, the distribution of jump 
rates contributes most to the non-exponentiality when 
the lithium ions successfully escape from their sites. 

The lifetime of the dynamical heterogeneities can be 
measured based on the four-time correlation function 



S 4 (ti2=t s ,t,t 34: = t s ) = <4(*12)4(*34.)> 



(9) 



In this experiment, a dynamical filter selecting lithium 
ions that are slow on the timescale is applied dur- 
ing two time intervals t\2 = t s and £34 = t s separated by 
a time t. Thus, S4(t Sl t,t s ) is given by the fraction of 
lithium ions that are slow during t 12 and a time t later 
during £34. Consequently, the signal decreases when slow 
ions become fast in the time interval t and mapping S4 (t) 
allows us to measure the timescale of the exchange of the 
dynamical state. In Fig. ^3 various normalized corre- 
lation functions S£{t) are compiled. In panel (a), we 
show results for different filter times t s at T = 713K 
and, in panel (b), we display data for comparable t s 
at various temperatures (S2{t s ) ~ 0.48). For the nor- 
malization, 5*4(7:) was fitted to a modified KWW func- 
tion, A 54 exp[— (t/r)P] + S4 , and scaled according to 
Sf(t) = (S 4 (t)- S^)/AS 4 . Evidently, S£{t) depends 
on both the filter time and the temperature where all 
decays are non-exponential. The non-exponentiality in- 
dicates that the exchange processes between slow and fast 
ions cannot be described by a single time constant, but 
are rather governed by a broad distribution of exchange 
times. Further, since dynamically diverse subensembles 
are selected for different filter times t s , the dependence on 
this parameter implies that the exchange time depends 
on the dynamical state itself. A comparison with Fig. [5] 
reveals that the decays of S2 (t) and S4 (t) are complete 
at comparable times. In accordance with the interpreta- 
tion of G4, this finding shows that the information about 
the initial dynamical states is lost when all lithium ions 
have successfully escaped from their sites. The mean 
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time constants T54, i.e., the mean exchange times, can 
be calculated from the parameters of the KWW fits to 
Si(t). Revisiting Fig. we see that ts4 and TS2 show 
a similar temperature dependence at T <T g . All these 
findings consistently indicate that the exchange of the dy- 
namical state in the glass results from successful jumps 
among neighboring lithium sites. Above T g , a stronger 
temperature dependence of tsa is observed. Hence, the 
modification of the lithium sites due to matrix dynamics 
provides another channel for the exchange of the dynam- 
ical state in the melt. 



IV. DISCUSSION 

In this article, we presented MD simulations performed 
to study the dynamics of non-crystalline LiPC>3 in a tem- 
perature range where the motions of the various atomic 
species strongly decouple, i.e., the dynamics in the melt 
(T >T g ) and in the glass (T <T g ) were compared. We 
found that the mechanism of the lithium motion changes 
from liquid-like diffusion to hopping dynamics when the 
temperature is decreased. Our main goal was to study 
the lithium ionic jumps at T <T g . Non-exponential two- 
time correlation functions, such as the incoherent inter- 
mediate scattering function, and sublinear diffusion over 
several orders of magnitude in time reflect the complexity 
of the lithium motion in a rigid phosphate glass matrix. 
For a comprehensive analysis of this dynamical process, 
we calculated suitable multi-time correlation functions. 
Three time-correlation functions showed that the non- 
exponential relaxation of the lithium ions at T <T g re- 
sults from both correlated back-and-forth jumps and the 
existence of dynamical heterogeneities, i.e., the presence 
of a distribution of jump rates. Measuring the lifetime 
of the dynamical heterogeneities by means of four-time 
correlation functions, we found that the non-uniformities 
of the motion are short-lived. In what follows, we sum- 
marize and discuss our results in more detail where the 
quantities that lead to the respective conclusion are in- 
dicated in brackets. 

The dynamical heterogeneities at T<T g show the fol- 
lowing features: (i) The distribution of jump rates con- 
tributes more to the non-exponentiality the lower the 
temperature (S3, H). (ii) For a given temperature, its 
contribution is maximum when the lithium ions success- 
fully escape from their sites (th/ts2)- (hi) The exchange 
between the fast and the slow ions of the distribution 
takes place on the same timescale as the lithium jumps, 
i.e., the dynamical heterogeneities are short-lived (G4, 
S4). In particular, both processes show a similar tem- 
perature dependence (jsa/ts2)- (iv) The timescale of 
the exchange depends on the dynamical state itself (S4). 
(v) The dynamical heterogeneities mainly result from dif- 
ferent jump rates at neighboring sites, but there are also 
non-uniformities of the motion on a length scale r > 10A 
(a x . y ). The back-and-forth jumps of the lithium ions in 
the glass can be characterized as follows: (i) Correlated 



back-and-forth jumps occur for the fast ions of the dis- 
tribution, but not for the slow ions, i.e., the back-jump 
probability depends on the initial dynamical state (G3, 
2^23)- (ii) For highly mobile ions, there is a backward cor- 
relation even beyond the first neighbor shell (^23). Well 
above T g , dynamical heterogeneities and back-and-forth 
motions on the length scale tmli are absent (5?23j <^x,y)- 
However, due to the temporary caging of the lithium ions 
by the neighboring particles, anisotropic dynamical het- 
erogeneities and a backward correlation exist on short 
length- and timescales. 

We conclude that the formation of well-defined lithium 
sites is accompanied by the emergence of correlated 
back-and-forth jumps and dynamical heterogeneities on 
a length scale r > run . At the lower temperatures, the 
non-exponential depopulation of the lithium sites results 
to a large extent from the dynamical heterogeneities and, 
hence, the distribution of jump rates is broad. The short 
lifetime of the dynamical heterogeneities implies that 
sites featuring fast and slow lithium dynamics, respec- 
tively, are intimately mixed. In particular, the absence of 
long-lived heterogeneities excludes a micro-segregation of 
the glass into extended clusters in which the lithium ions 
are mobile and immobile, respectively. Instead, we sug- 
gest that the broad distribution of jump rates in LiPC>3 
glass results from a diversity of the local glass structure 
at neighboring sites, e.g., from varying coordination num- 
bers. In this basically random energy landscape, fast 
ion migration occurs along pathways connecting low en- 
ergy barriers. The presence of such channels is not only 
indicated by the appearance of the lithium trajectories, 
but also by the existence of long-range dynamical het- 
erogeneities and by the backward correlation beyond the 
first neighbor shell for the fast ions. 

Some of our findings are recovered when analyzing 
the dynamics resulting from schematic models such as 
the random-barrier and the random-energy model where 
the particles move in a time-independent, external en- 
ergy landscape. For example, the particles migrate along 
preferential pathways in these models and percolation 
approaches have been applied to describe the dynam- 
ics at low temperatures£ii£2iS. Moreover, the back-jump 
probability depends on the dynamical state of a particle, 
e.g., for the random-barrier model, many back-and-forth 
jumps occur over low energy barriers, whereas slow ions 
cross high barriers so that a backward motion is unlikely. 
However, other features of the lithium dynamics in LiPC>3 
glass are not reproduced in schematic models. At vari- 
ance with our findings, cf. Fig. 1121 a straightforward cal- 
culation shows that H it — ► 0) ~ 0.5 is obtained for these 
models. Furthermore, we observe that the first maxi- 
mum of o x ^ v {t\2) lies at r > ruu/2 and shifts towards 
shorter distances when t is increased, cf. Fig. This 
means that, in contrast to what is captured in simple 
models, the energy barriers crossed by the fast lithium 
ions of the distribution are located at larger distances 
than the ones overcome by the slow ions. We speculate 
that simple hopping models are suited to study the time- 
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independent effects of the glassy matrix on the lithium 
ionic diffusion, but they fail in modelling the influence 
of the time-dependent Li-Li interactions. For example, 
we want to demonstrate that our findings for the posi- 
tion of the energy barrier can be explained when assum- 
ing that the Li-Li interaction decreases with the waiting 
t w at a site due to the adjustment of the neighboring 
lithium ions. For simplicity, we consider a lithium ion 
that occupies a site at r = and experiences a poten- 
tial V(r,t w ) = —cos(nr) + (l/2)a(t w )r 2 . Here, the first 
and the second term are meant to mimic the contribu- 
tion from the glassy matrix and from the repulsion due 
to other lithium ions (a(t w ) > 0), respectively. A Tay- 
lor expansion for r rs 1 shows that the maximum of this 
potential lies at r m (t w ) = 1/(1 — a(t w )/n 2 ). Hence, for 
da(t w ) I dt w < 0, this simple model yields a maximum at 
r m (tw) > 1 that is located at smaller distances the larger 
the value of the waiting time as observed in our simula- 
tions. 

Comparing the results of MD simulation studies on 
phosphate and silicate glasses several similar features 
of cation dynamics are striking: (i) It is widely ac- 
cepted that cation diffusion at sufficiently low tem- 
peratures can be described as a sequence of hopping 
processea 19 ! 35 ! 39 ! 41 ! 42 ! 43 ! 44 . (ii) F s (q m ,t) consistently 
shows non-exponential decays that are well described by 
KWW functions with stretching parameters (3 — 0.42 — 
0.47^i2i. (iii) Dynamical heterogeneities with a limited 
lifetime exists that become more pronounced when the 
temperature is decreased 1 ^, (iv) Fast ion transport along 
preferential pathways takes plac o 22 ' 24 . (v) A high back- 
jump probability exists for the fast ions of the distribu- 
tion, but not for the slow ions 4 ^. On the other hand, 
in previous work on sodium silicate glasses^, correlated 
back-and-forth jumps were not observed at somewhat 
higher temperatures. However, since, a high back-jump 
probability exists only for the fast ions at sufficiently low 
temperatures, the percentage of ions showing correlated 
back-and-forth jumps may have been too small to be de- 
tected in the study by Sunyer et alM>. Altogether, though 
it is to be further clarified to which extent the findings of 
MD simulations reflect the dynamics in a specific system 
at lower temperatures, the similarity of the simulation re- 
sults for various interaction potentials suggests that the 
mechanism of cation dynamics in various glassy ion con- 
ductors is highly comparable. 

Finally, we compare our results with experimental find- 
ings. First, it can be noted that two-time correlation 
functions of lithium dynamics obtained in experiments 
and simulations for LiP03 glass show a similar non- 
exponentiality. Specifically, the stretching parameter /3 ~ 
0.42 found here compares well to = 0.33—0.54 observed 
in electrical and mechanical relaxation studies 3 ' 4 5 ' 12 and 
in multi-dimensional 7 Li NMR experiments? 4 . Further, 
in these measurements, an activation energy E a «0.68eV 
was obtained for the lithium dynamics. Comparing this 
value with the results of our constant pressure simula- 
tions good agreement is found for Eps ~0.67eV, but sig- 



nificant differences must be noted for Ep, = 0.47eV. For 
the moment, it is not clear where the difference between 
Eps and Er> results from. Since similar deviations were 
reported in MD simulation studies of LiSi02 glass 1 ^, one 
may speculate that they are a consequence of the com- 
plexity of ion dynamics in glasses. 

In addition, our findings are useful for the interpre- 
tation of experimental results on the origin of the non- 
exponential relaxation of the mobile ions in solid ion con- 
ductors. Multi-dimensional 109 Ag NMR experiments on 
glassy and crystalline silver ion conductors'^ 5 - showed 
that the non-exponential silver ionic relaxation is due to 
the existence of a broad distribution of jump rates rather 
than to correlated back-and-forth jumps. On the other 
hand, the presence of correlated back-and-forth motion 
is indicated by a frequency dependence of the electri- 
cal conductivity q^WWW 7 . We suggest that 
this apparent discrepancy can be cleared up based on 
the present results. Specifically, we found that, first, the 
back-jump probability depends on the dynamical state, 
see alsoi 5 -, and, second, the contribution of the dynam- 
ical heterogeneities to the non-exponentiality increases 
upon cooling. Since cr(cj) is dominated by the fast ions, 
we expect a disproportional contribution from the cor- 
related back-and-forth jumps of ions in this dynamical 
state. In contrast, multi-dimensional NMR experiments 
probe the behavior of non-fast ions at low temperatures 
and, hence, if at all, back-and-forth jumps should be of 
minor relevance. 

In summary, various techniques have to be applied si- 
multaneously to obtain a complete picture of a complex 
dynamical process such as ion dynamics in solids. In par- 
ticular, a further comparison of results obtained in MD 
simulations and in multi-dimensional NMR experiments 
is very promising since both techniques allow one to study 
a dynamical process on a microscopic level. For such a 
comparison, the multi-time correlation functions S2, S3 
and S4 can be particularly useful, because their informa- 
tion content is similar to the one of the corresponding 
quantities in NMR experiments. 
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APPENDIX A: CALCULATION OF 
THREE-TIME CORRELATION FUNCTIONS 

It is well established for supercooled liquids that, for 
purely homogeneous and purely heterogeneous relax- 
ation, three-time correlation functions comparable to S3 
can be expressed by the corresponding two-time corre- 
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lation functiono 18 ! 46 ! 48 ! 49 . Here, S 3 (t/2, t/2) = [S 2 {t/2)] 2 
follows, because no selection is possible for purely homo- 
geneous dynamics. Hence, the dynamics during ti 2 and 
£23, respectively, are statistically independent so that the 
joint probability for an escape in neither time interval 
can be written as a product. If this was not true a dy- 
namically diverse subensemble could be selected based 
on the motion in the first time interval and, thus, dy- 
namical heterogeneities would exist. For purely hetero- 
geneous dynamics, the correlation functions result from 
a superposition of the contributions of exponentially re- 
laxing subensembles. For the subensemble characterized 
by a correlation time r^, one has 



S 2 ,i(t) = exp - 



t/2 + 1/2 



= S 2 4t/2)S 2 4t/2) 



= S 3)i (t/2,t/2). 



The second equals sign is due to the definition of the 
exponential function, and the third results, because the 
relation for purely homogeneous dynamics holds for the 
exponentially relaxing, undecomposable subensembles. 
Since S 3 ^{t/2,t/2)=S 2 ^{t) is valid for all subensembles, 
S 3 (t/2,t/2) = S 2 (t) follows. 
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